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Abstract 

Effects of the propagation of particles, which have a finite hfe-time and an ac- 
cording width in their mass spectrum, are discussed in the context of transport de- 
scription. First, the importance of coherence effects (Landau-Pomeranchuk-Migdal 
effect) on production and absorption of field quanta in non-equilibrium dense matter 
is considered. It is shown that classical diffusion and Langevin results correspond 
to re-summation of certain field-theory diagrams formulated in terms of full non- 
equilibrium Green's functions. Then the general properties of broad resonances in 
dense and hot systems are discussed in the framework of a self-consistent and con- 
serving <I>-derivable method of Baym at the examples of the p-meson in hadronic 
matter and the pion in dilute nuclear matter. Further we address the problem 
of a transport description that properly accounts for the damping width of the 
particles. The ^-derivable method generalized to the real-time contour provides a 
self-consistent and conserving kinetic scheme. We derive a generalized expression 
for the non-equilibrium kinetic entropy fiow, which includes corrections from fluc- 
tuations and mass-width effects. In special cases an ff-theorem is proved. Memory 
effects in collision terms give contributions to the kinetic entropy flow that in the 
Fermi-liquid case recover the famous bosonic type In T correction to the specific 
heat of liquid Helium-S. At the example of the pion-condensate phase transition in 
dense nuclear matter we demonstrate important part played by the width effects 
within the quantum transport. 

Quasiparticle representations in many-body theory were designed by Landau, Migdal, 
Galitsky and others, see refs [0, 0, This concept was first elaborated at the example 
of low-lying particle-hole excitations in Fermi liquids. A.B. Migdal was the first to apply 
these methods to description of various nuclear phenomena and construction of a closed 
semi-microscopic approach that is now usually called "Theory of finite Fermi systems" 
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The need for explicit treatment of soft modes within this approach stimulated him to 
generalize this concept to include pion and A excitations. A.B. Migdal suggested variety 
of interesting effects like softening of the pion mode in nuclei, pion condensation in dense 
nuclear and neutron star matter and possible existence of abnormal nuclei glued by pion 
condensate 0, ||. These ideas stimulated further development of pion physics 

with applications to many phenomena in atomic nuclei, neutron stars and heavy-ion 
collisions, see |l^ and refs therein. In this paper we would like to briefly review 

recent developments of some of the above ideas as applied to heavy-ion physics. 

With the aim to describe the collision of two nuclei at intermediate and high energies 
one is confronted with the fact that the dynamics has to include particles like the in- 
medium excitation with the pion quantum numbers, as well as the delta and rho-meson 
resonances with life-times of less than 2 fm/c or equivalently with damping widths above 
100 MeV. Also the collisional damping rates deduced from presently used transport codes 
are of the same order, whereas typical mean temperatures range between 50 to 150 MeV 
depending on beam energy. Thus, the damping width of most of the constituents in 
the system can by no means be treated as a perturbation. As a consequence the mass 
spectrum of the particles in the dense matter is no longer a sharp delta function but 
rather acquires a width due to collisions and decays. One thus comes to a picture which 
unifies resonances, which have already a width in vacuum due to decay modes, with the 
"states" of particles in dense matter, which obtain a width due to collisions (collisional 
broadening). Landau, Pomeranchuk and Migdal were first to demonstrate importance of 
finite-width effects in multi-particle scattering |Tl], |T^. Such a coherence scattering effect, 
known now as Landau-Pomeranchuk-Migdal effect, was recently observed at Stanford 



accelerator 13 



The theoretical concepts for a proper many-body description in terms of a real time 
non-equilibrium field theory have been devised by Schwinger [|T^ , Kadanoff and Baym |T5 



and Keldysh |jT6| in the early sixties. However a proper dynamical scheme in terms of 
a transport concept that deals with unstable particles, like resonances, is still lacking. 
Rather ad-hoc recipes are in use that sometimes violate basic requirements as given by 
fundamental symmetries and conservation laws, detailed balance or thermodynamic con- 
sistency. The problem of conserving approximation has first been addressed by Baym and 
Kadanoff [|1^, [1^ . They started from an equilibrium description in the imaginary-time for- 
malism and discussed the response to external disturbances. Baym, in particular, showed 
Tsj that any approximation, in order to be conserving, must be so-called ^-derivable. It 



turned out that the $ functional required is precisely the auxiliary functional introduced 
by Luttinger and Ward ||l9l (see also ref. ||2^) in connection with the thermodynamic 
potential. In non-equilibrium case the problem of conserving approximations is even more 
severe than in near-equilibrium linear-response theory |2l|, |2^ . 

In this review we discuss recent developments of the transport theory beyond the 
quasiparticle approximation and consequences of the propagation of particles with short 
life-times in hadron matter. First we consider few examples of equilibrium systems which 
clearly indicate that treatment beyond the quasiparticle approximation is really needed. 
We start with a genuine problem related to the occurrence of broad damping width, i.e. 
the soft-mode problem (Landau-Pomeranchuk-Migdal effect). This is the direct radiation 
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of quanta from a piece of a dense medium [Q. Classically this problem is formulated 
as coupling of a coherent classical field, e.g., the Maxwell field, to a stochastic source 
described by the Brownian motion of a charged particle. In this case the classical current- 
current correlation function, can be obtained in closed analytical terms and discussed as 
a function of the macroscopic transport properties, the friction and diffusion coefficient 
of the Brownian particle. We show that this result corresponds to a partial re-summation 
of photon self-energy diagrams in the real-time formulation of field theory. Subsequently, 
properties of particles with broad damping width are illustrated at the examples of the 
p-meson in dense matter at finite temperature ||2^ and the pion in the limit of a dilute 



nuclear matter . The question of consistency becomes especially important for a multi- 
component system like nNAp, where the properties of one species can change due to 
the presence of interactions with the others and vice versa. The "vice versa" is very 
important and corresponds to the principle of actio = re-actio. This implies that the 
self-energy of one species cannot be changed through the interaction with other species 
without affecting the self-energies of the latter ones also. The <l>-derivable method of 
Baym []TB| offers a natural and consistent way to account for this principle. 

Then we address the question how particles with broad mass-width can be described 



consistently within a transport picture [21, 22|. We argue that the Kadanoff-Baym 



equations in the first gradient approximation together with the ^-functional method of 
Baym |jl8| provide a proper self-consistent approach for kinetic description of systems of 
particles with a broad mass-width. We argue that after gradient expansion the full set 
of equations describing quantum transport contains two equations, the differential gen- 
eralized kinetic equation for a distribution function in 8-(X, p)-space and the algebraic 
equation for the spectral density. Other equations are resolved. We discuss the problems 
of concerning charge and energy-momentum conservation, thermodynamic consistency, 
memory effects in the collision term and the growth of entropy in specific cases. Finally, 
we demonstrate finite-width effects in quantum kinetic description at the example of pion 
condensation, where the width of soft pionic excitations due to their decay into particle- 
hole pairs governs the dynamics of the phase transition in the isospin-symmetric nuclear 
matter. 

We use the units Ti = c = 1. For simplicity we treat fermions non-relativistically 
whereas bosons (mesons) are treated as relativistic particles. 



1 Bremsstrahlung from Classical Sources 



For a clarification of the soft mode problem following |2^ we first discuss an example in 
classical electrodynamics. We consider a stochastic source, the hard matter, where the 
motion of a single charge is described by a diffusion process in terms of a Fokker-Planck 
equation for the probability distribution / of position x and velocity v 

Fluctuations also evolve in time according to this equation, or equivalently by a random 
walk process and this way determine correlations. This charge is coupled to the 
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Figure 1: Left: Current-current correlation function in units of < f ^ > as a function of 
time (in units of l/F^.) for different values of the photon momentum = 3A;^r^/< > 
with A; = 0, 1, 2, 3. Right: Rate of real photons d^N/{dujdt) in units of 47re^ {v"^) /3 for a 
non-relativistic source for =50,100,150 MeV; for comparison the IQF results (dashed 
lines) are also shown. 



Maxwell field. On the assumption of a non-relativistic source, this case does not suffer 
from standard pathologies encountered in hard thermal loop (HTL) problems of QCD, 
namely the coUinear singularities, where vq ^ 1, and from diverging Bose-factors. The 
advantage of this Abelian example is that damping can be fully included without violat- 
ing current conservation and gauge invariance. This problem is related to the Landau- 



Pomeranchuk-Migdal effect of bremsstrahlung in high-energy scattering [0, [12 



The two macroscopic parameters, the spatial diffusion D and friction coefficients 
determine the relaxation rates of velocities. In the equilibrium limit (t oo) the 
distribution attains a Maxwell-Boltzmann velocity distribution with the temperature 
T = m {v"^) /3 = mDTx- The correlation function can be obtained in closed form and one 
can discuss the resulting time correlations of the current at various values of the spatial 
photon momentum q. Fig. |l| (details are given in ref. [23|). For the transverse part of the 



correlation tensor this correlation decays exponentially as ~ e~^'^'^ at q = 0, and its width 
further decreases with increasing momentum q = \q\. The in-medium production rate is 
given by the time Fourier transform r — > cj. Fig. |I| (right part). The hard part of the 
spectrum behaves as intuitively expected, namely, it is proportional to the microscopic 
collision rate expressed through (cf. below) and thus can be treated perturbatively by 
incoherent quasi-free (IQF) scattering prescriptions. However, independently of F^. the 
rate saturates at a value of ~ 1/2 in these units around uj ~ Fj,, and the soft part shows 
the inverse behavior. That is, with increasing collision rate the production rate is more 
and more suppressed! This is in line with the picture, where photons cannot resolve the 
individual collisions any more. Since the soft part of the spectrum behaves like cu/Fa,, it 
shows a genuine non-perturbative feature which cannot be obtained by any power series 
in Fa;. For comparison: the dashed lines show the corresponding IQF yields, which agree 
with the correct rates for the hard part while completely fail and diverge towards the soft 
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Figure 2: Self-energy diagrams determining the production and absorption rates. 



end of the spectrum. For non-relativistic sources (v^) <C 1 one can ignore the additional 
g-dependence (dipole approximation; cf. Fig. |l|) and the entire spectrum is determined 
by one macroscopic scale, the relaxation rate F^,. This scale provides a quenching factor 



Coiu) = ^— - (2) 



Y2 



X 



by which the IQF results have to be corrected in order to account for the finite collision 
time effects in dense matter. 

The diffusion result represents a re-summation of the microscopic Langevin multiple 
collision picture and altogether only macroscopic scales are relevant for the form of the 
spectrum and not the details of the microscopic collisions. Note also that the classical 



result fulfills the classical version (h 0) of the sum rules discussed in refs. |23 



2 Radiation at the Quantum Level 

We have seen that at the classical level the problem of radiation from dense matter 
can be solved quite naturally and completely at least for simple examples, and Figs. ^ 
display the main physics. They show, that the damping of the particles due to scattering 
is an important feature, which in particular has to be included right from the onset. 
This does not only assure results that no longer diverge, but also provides a systematic 
and convergent scheme. On the quantum level such problems require techniques beyond 
the standard repertoire of perturbation theory or the quasi-particle approximation. In 
terms of non-equilibrium diagrammatic technique in Keldysh notation, the production 
or absorption rates are given by the self-energy diagrams of the type presented in Fig. 
2 with an in- and out-going radiating particle (e.g. photon) line The hatched 

loop area denotes all strong interactions of the source. The latter give rise to a whole 
series of diagrams. As mentioned, for the particles of the source, e.g. the nucleons, one 
has to re-sum Dyson's equation with the corresponding full complex self-energy in order 
to determine the full Green's functions in dense matter. Once one has these Green's 
functions together with the interaction vertices at hand, one could in principle calculate 
the required diagrams. However, both the computational effort to calculate a single 
diagram and the number of diagrams increase dramatically with the loop order of the 
diagrams, such that in practice only lowest-order loop diagrams can be considered in the 
quantum case. In certain limits some diagrams drop out. We could show that in the 
classical limit, which in this case implies the hierarchy a;, |qr|, F ^ T ^ m together with 
low phase-space occupations for the source, i.e. f{x,p) <^ 1, only the following set of 



diagrams survives 



+ 



+ + 



+ 



(3) 



In these diagrams the bold hnes denote the full nucleon Green's functions which also 
include the damping width, the black blocks represent the effective nucleon-nucleon inter- 
action in matter, and the full dots, the coupling vertex to the photon. The thin dashed 
lines show how the diagrams are to be cut into two interfering amplitudes. This way each 
of these diagrams with n interaction loop insertions just corresponds to the n*^ term in the 
corresponding classical Langevin process, where hard scatterings occur at random with a 
constant mean collision rate T . These scatterings consecutively change the velocity of a 
point charge from to Wi, to 1)2, • • •• In between scatterings the charge moves freely. For 
such a multiple collision process the space integrated current-current correlation function 
takes a simple Poisson form 



oc j d3sid3a;2(j"(a;i, t - |)j'^(a;2, t + 



E 

n=0 



(4) 



with V = Here (. . .) denotes the average over the discrete collision sequence. This 

form, which one writes down intuitively, agrees with the analytic result of the quantum 
correlation diagrams (0) in the limit f{x,p) <^ 1 and T <^ T. Fourier transformed it 
determines the spectrum in completely regular terms (void of any infra-red singularities), 
where each term describes the interference of the photon being emitted at a certain time or 
n collisions later. In special cases where velocity fluctuations are degraded by a constant 
fraction a in each collision, such that {vq ■ Vn) = (f ■ ^o); one can re-sum the whole 
series in Eq. (1) and thus recover the relaxation result with (v'^) = F {{vq — viY) at 
least for g = and the corresponding quenching factor (^. Thus the classical multiple 
collision example provides a quite intuitive picture about such diagrams. Further details 
can be found in ref. . 

The above example shows that we have to deal with particle transport that explicitly 
takes account of the particle mass-width in order to properly describe soft radiation from 
the system. 



3 The p-meson in Dense Matter 

Following $ derivable scheme we will first discuss two examples within thermo-equilibrium 
systems. First we will concern properties of the p-meson and their consequences for the 



p-decay into di-leptons In terms of the non-equilibrium diagrammatic technique, the 
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exact production rate of di-leptons is given by the following formula 



^j^e+e 
dtdm 




= fp{m,p,x,t) Ap{m,p,x,t) 2mTf"^^'^ (m). (5) 

Here r^^^^ (m) oc is the mass-dependent electromagnetic decay rate of the p into 

the di-lepton pair of invariant mass m. The phase-space distribution fp{m, p, x, t) and the 
spectral function Ap{m,p,x,t) define the properties of the p-meson at space-time point 
X, t. Both quantities are in principle to be determined dynamically by an appropriate 
transport model. However till to-date the spectral functions are not treated dynamically 
in most of the present transport models. Rather one employs on-shell (5-functions for all 
stable particles and spectral functions fixed to the vacuum shape for resonances. 

As an illustration, we follow the two channel example presented by one of us in ref. . 



There the p-meson just strongly couples to two channels, i.e. the vr+vr" and nN ^ pN 
channels, the latter being relevant at finite nuclear densities. The latter component is 
representative for all channels contributing to the so-called direct p in transport codes. For 
a first orientation the equilibrium properties^ are discussed in simple analytical terms with 
the aim to discuss the consequences for the implementation of such resonance processes 
into dynamical transport simulation codes. 

Both considered processes add to the total width of the p-meson 

Ttot{m,p) = Tp^^+^-{m,p) + Tp^^NN-'^{m,p), (6) 

and the equilibrium spectral function then results from the cuts of the two diagrams 

TT , N 




2mrp^+^- +2mrp^jvjv^i 
(m2 - m2 - ReS^)^ + rv?Yl^^ 

In principle, both diagrams have to be calculated in terms of fully self-consistent propa- 
gators, i.e. with corresponding widths for all particles involved. This formidable task has 
not been done yet. Using micro- reversibility and the properties of thermal distributions. 



^Far more sophisticated calculations have already been presented in the literature ^ ^ 
]. It is not the point to compete with them at this place. 
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the two terms in Eq. contributing to the di-lepton yield (^), can indeed approximately 
be reformulated as the thermal average of a tt+tt^ ^ p — e"'"e~-annihilation process and 
a T[N — > pN e"'"e~A^-scattering process, i.e. 

oc (U+U- v^^ o-(7r+7r" ^ e+e") 



dmdt 

+ fnfN v^N (r{7iN ^pN^ e+e-iV))^ , (8) 

where f-,, and are corresponding particle occupations and Vt,t, and Vt^n are relative 
velocities. However, the important fact to be noticed is that in order to preserve unitarity 
the corresponding cross sections are no longer the free ones, as given by the vacuum 
decay width in the denominator, but rather involve the medium dependent total width 
(^) . This illustrates in simple terms that rates of broad resonances can no longer simply 
be added in a perturbative way. Since it concerns a coupled channel problem, there 
is a cross talk between the different channels to the extent that the common resonance 
propagator attains the total width arising from all partial widths feeding and depopulating 
the resonance. While a perturbative treatment with free cross sections in Eq. @ would 
enhance the yield at resonance mass, m = rrip, if a channel is added, cf. Fig. 3 left 
part, the correct treatment (0) even inverts the trend and indeed depletes the yield at 
resonance mass, right part in Fig. 3. Furthermore, one sees that only the total yield 
involves the spectral function, while any partial cross section refers to that partial term 
with the corresponding partial width in the numerator! Compared to the shape of the 
spectral function both thermal components in Fig. 3 show a significant enhancement 
on the low mass side and a strong depletion at high masses due to the thermal weight 
/ oc exp(— po/r) in the rate (j^). This kinematical effect related to the broad width also 
survives in non-equilibrium calculations and is a signature of phase space restrictions 
imposed for particles with higher energies. The related question how to preserve detailed 
balance in the case of broad resonances was addressed by Danielewicz and Bertsch . 
The latter method has then been implemented in transport models mostly applied to 
the description of the A-resonance. For the transport description of the p meson only 
quite recently a description level has been realized that properly includes the width effects 



discussed above, e.g. in ref. |3g], cf. also the comments given in ||33]. The transport 
treatment of broad resonances is further discussed in sections 4-7. 

As an example we show an exploratory study of the interacting system of vr, p and 
Oi-mesons described by the ^-functional 




$ = i-N^N^N^N^sA + ^.✓N^N^N^V* + 7]- (9) 



(cf. section ^ below). The couplings and masses are chosen as to reproduce the known 
vacuum properties of the p and ai meson with nominal masses and widths rrip = 770 MeV, 
TUa, = 1200 MeV, Tp = 150 MeV, F^, = 400 MeV. The results of a finite temperature 
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Di-lepton rates from thermal p-mesons (T = 110 MeV) 




Figure 3: e"'"e~ rates (arb. units) as a function of the invariant pair mass m at T = 110 
MeV from tt'^tt'^ annihilation (dotted line) and direct p-meson contribution (dashed line), 
the full line gives the sum of both contributions. Left part: using the free cross section 
recipe, i.e. with Ftot = ^p-k+tt-', right part: the correct partial rates (0). Ap is given by the 
thick line. The calculations are done with Tp^.„.„{mp) = 150 MeV and Tp^T^]\rj^-i{mp) = 70 
MeV. 



Dilepton rate (arb.units), T=110 MeV, p=150 MeV/c 




inv. Mass (MeV) 

Figure 4: Left part: contour plot of the self-consistent spectral function of the p-meson 
as a function of energy and spatial momentum; right part: thermal di-lepton rate as a 
function of invariant mass at | p | = 300 MeV/c 
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calculation at T = 150 MeV with all self-energy loops resulting from the ^-functional of 
Eq. (H) computed with self-consistent broad width Green's functions are displayed in 
Fig. 4 (corrections to the real part of the self-energies were not yet included). The last 
diagram of $ with the four pion self-coupling has been added in order to supply pion with 
broad mass-width as they would result from the coupling of pions to nucleons and the 
A resonance in nuclear matter environment. As compared to first-order one-loop results 
which drop to zero below the 2-pion threshold at 280 MeV, the self-consistent results 
essentially add in strength at the low-mass side of the di-lepton spectrum. 



3.1 Virial Limit 

In the dilute-density limit (virial limit) the corresponding self-energies of the particles 
and intermediate resonances are entirely determined by two-body scattering properties, in 
particular, by scattering phase shifts. We illustrate this at the example of the interacting 
system of nucleons, pions and delta resonances, which has recently been investigated by 
Weinhold et al. Following their study we consider a pedagogical example, where 

the TriViV-interaction is omitted. Then with a p-wave TriVA-coupling vertex among the 
three fields the first and only diagram of $ up to two vertices and the corresponding three 
self-energies are given by 





(10) 



Here the solid, dashed and double lines denote the propagators of A^, vr and A, respectively. 
In non-relativistic approximation for the baryons we ignore contributions from the baryon 
Dirac-sea. Then the bare pion mass agrees with its vacuum value, while the nucleon and 
delta masses require appropriate mass counter terms. The A self-energy Sa attains the 
vacuum width and position of the delta resonance due to the decay into a pion and a 
nucleon. The corresponding scattering diagrams are obtained by opening two propagator 
lines of $ with the prominent feature that the vriV-scattering proceeds through the delta 
resonance. Since in this case a single resonance couples to a single scattering channel, the 
vacuum spectral function of the resonance can directly be expressed through the scattering 
T-matrix and hence through measured scattering phase shifts 



33 



4sin^ (p) 



111 



where p = Pn + P-k- Thus through (0) the vacuum properties of the delta can almost 
model-independently be obtained from scattering data. 

For the multi-component system the renormalized thermodynamic potential including 
vacuum counter terms can be written as 



Vt {G^, Gn, Ga} 
= T E 

ae{7r,Af,A} 



T«:Tr {- In [-Gf (po + iO,p)] + GfSf }^^^ + $r 



(12) 
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1/2 for neutral bosons, 
K={ ' (13) 
1 for charged bosons and fermions. 

For any function f{p) the thermodynamic trace Tr{...}T^^ is defined as 

T^{/(p)}t,m = Td"^ I -0j-APo - P) 2 Im/(po + iO,p) (14) 

an energy integral over thermal occupations n(e) = [exp(£:/T) ± 1]"^, of Fermi-Dirac/Bose- 
Einstein type. The upper sign appears for fermions the lower, for bosons, d is the degen- 
eracy in that particle channel, and V denotes the volume. Eq. ([l^) still has the functional 
property to provide the retarded Dyson's equations for the from the stationary con- 
dition which we use in order to determine the physical value of Q. For the particular case 
here one further can exploit the property 

$r = ±nT TT{^aGa}T,^,■, for a e {N, n, A} and $ of form (0), 

valid for that depend linearly on all propagators. Compatible with the low density 
limit one can expand the Trln{— G} terms for the pion and nucleon around the free 
propagators, and finally obtains 



= + + T Tr {in [-G^(po + iO,p)] }^ (15) 



= + + d^TV j (^2^1^ In [1 - (po - /ia)] (16) 

for the physical value of Q. Here the ^1^^^° are the free single-particle thermodynamic 
potentials^, while fJ^A = Pn + Pn and d^ = 16 are the chemical potential and degeneracy 
factor of the A resonance, respectively. The last term in (|TB|) obtained through (|TT|), rep- 
resents a famous result derived by Beth-Uhlenbeck ||38| , |39| , later generalized by Dashen, 
Ma and Bernstein |Q and applied to nuclear resonance matter in refs. [0, ^ It 



illustrates that the virial corrections of the system's level density due to interactions are 
entirely given by the energy variation of the corresponding two-body scattering phase 
shifts d6/dpo. 

All thermodynamic properties can be obtained from Q through partial differentiations 
with respect to T and the /i. The final form (|T6]) may give the impression that one deals 
with non-interacting nucleons and pions. This is however not the case. For instance the 
densities of baryons and pions derived from (|I6|) become^ 

d^nNA free , , f^^yrAf A free , , / 1 ^7^ 

PB = ^ = Pn +PA+Pcorr, Ptt = ^ =P^ + PA + Pcorr, (17) 

OPN OPn 



^The appropriate cancellation of terms for the result jisl ) is only achieved, if one uses fl^'^'^'^, i.e. the 
partition sum of free particles with the free energy-momentum dispersion relation. Within this model 
already on the vacuum level the nucleon would acquire loop corrections to its self-energy which would 
lead to deviations between fi™'^ and ri^''®", as well as between the corresponding propagators off their 
mass shell. 

•^In equilibrium /i^r has to be put to zero after differentiation. 
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with 



f d^p 

pA = dA J 7^;^^a(P0 - /iA)^r^(p), (1^ 



(27r) 

J^^nj^{pO - fl/^)Bcorr{p), (19) 

Here the density of dehas is determined by the delta spectral function. The interaction 
contribution contained in the correlation density pcorr depends on the difference between 
the phase-shift variation and the spectral function 



5_ = 2^%^^ - ATip) = 2Im 
dpo 



dpo 



(20) 



Due to the fact that Ta{p) grows with energy and the real part of Ga changes sign at the 
resonance energy, -Bcorr becomes positive below and negative above resonance, respectively. 
It leads to an enhancement of both densities at low energies, i.e. below resonance and 
this way to a further softening of the resulting equation of state compared to the naive 
spectral function treatment ignoring the -Bcorr terms. This illustrates that an interacting 
resonance gas cannot consistently be described by a set of free particles (here the pions and 
nucleons) plus vacuum resonances (here the delta), described by their spectral function. 
Rather the coupling of a bare resonance to the stable particles determines its width, and 
thus its spectral properties in vacuum. At the same time the stable particles are modified 
due to the interaction with the resonance. Only the account of all three self-energies in 
(p!0|) provides a conserving and thermodynamically consistent approximation. 

Alternatively to the picture above, the properties of the system can be discussed 
entirely in terms of the asymptotic particles, i.e. the pion and the nucleon. The thermo- 
dynamic potential is then still given by ([T6|). This form is valid even without intermediate 
resonances and the phase-shifts just account for the tcN interaction properties. Also the 
self-energy of the pion can be obtained from phase shifts by means of the optical theorem 



To linear order in the nucleon density pj\f one determines the pion self-energy 



S.(piab) = 4vr^p^F.^(0) = -27r^^ ^2 sinfee^'-, (21) 

Pcm Pcm ""ir 

from the forward vrA^'-scattering amplitude Ft,n{0). Here pi^h and pcm refer to the pion 
3-momenta in the matter rest frame and the cm. frame of the irN collisions. The arising 
kinematical factor piab/Pcm = y/s/^N which has mostly escaped notice even in standard 



references on the low desity theorem, e.g. becomes important for heavier projectiles 



like kaons, cf. ref. [^. Here the degeneracy factors dN '■ d.„ : = 4: : 3 : 16 just provide 
the proper spin/isospin counting. This self-energy, which determines an optical potential 
or index of refraction, is attractive below the delta resonance energy and repulsive above. 
It agrees with a related effect in optics, where a resonance in the medium causes an 
anomalous behavior of the real part of the index of refraction, which is larger than 1 
below the resonance frequency and less than 1 above the resonance. Thus, absorption, 
e.g. by exciting a resonance, is always accompanied by a change of the real part of the 
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index of refraction of the scattered particle. The ^-derivable principle automatically takes 
care about these features. 

As has been discussed in | J7|, the corrections to the system's level density (last term in 



(p!6|)) can also be inferred from the time shifts (or time delays) induced by the scattering 



processes. From ergodicity arguments one obtains for a single partial wave 

_d_ 

dpo 



d 

a (^lcvel(Po) ~ ^kvol(Po)) = ^forward + ^scatt. = ^delay 



= [sin 533 COS 533] +4 sin^ 533^ = 2^. (22) 
opo dpo dpo 

These expressions apply to the cm. frame. Here the forward delay time Tforward relates 
to the change of the group velocity induced by the real part of the optical potential, c.f. 



(pT|). The scattering time Tgcatt. finally results from the delayed re-emission of the pion 



from the intermediate resonance to angles off the forward direction. 



4 Quantum Kinetic Equation 

The three above-presented examples unambiguously show that for a consistent dynamical 
treatment of non-equilibrium evolution of soft radiation and broad resonances we need 
a transport theory that takes due account of mass-widths of constituent particles. A 
proper frame for such a transport is provided by Kadanoff-Baym equations. We consider 
the Kadanoff-Baym equations in the first-order gradient approximation, assuming that 
time-space evolution of a system is smooth enough to justify this approximation. 

First of all, it is helpful to avoid all the imaginary factors inherent in the standard 
Green's function formulation {G^^ with i,j G { — h}) and introduce quantities which are 
real and, in the quasi-homogeneous limit, positive and therefore have a straightforward 
physical interpretation [^], much like for the Boltzmann equation. In the Wigner repre- 
sentation we define 

F{X,p) = A{X,p)f{X,p) = (T)iG-+(X,p), 

F{X,p) = AiX,p)[l T /(X,p)] = iG+-(X,p), (23) 
A{X, p) = -21mG^{X, p)=F±F = i (^G+- - G'"+) (24) 

for the generalized Wigner functions F and F with the corresponding /onr-phase-space 
distribution functions f{X,p) and Fermi/Bose factors [1 =F f{^,p)], with the spectral 
function A{X,p) and the retarded propagator G^. Here and below the upper sign corre- 
sponds to fermions and the lower one, to bosons. According to relations between Green's 
functions G*-' only two independent real functions of all the G'-' are required for a com- 
plete description. Likewise the reduced gain and loss rates of the collision integral and 
the damping rate are defined as 

TUX,p) = (T)iS-+(X,p), T,^t{X,p) = iS+-(X,p), (25) 
T{X,p) = -2ImS^(X,p) = T,^t{X,p)±TUX,p), (26) 
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where S*-' are contour components of the self-energy, and is the retarded self-energy. 

In terms of this notation and within the first-order gradient approximation, the Kadanoff- 
Baym equations for F and F (which result from differences of the corresponding Dyson's 
equations with their adjoint ones) take the kinetic form 



VF - {lin, ReG^j = C, 
PF-{rout,ReG^} = TC^ 
with the drift operator and collision term respectively 



^ (9ReE^\ _ (9ReE^ d 
V = I 2n„ ^ I + 



C{X,p) = TUX,p)F{X,p) - T,^,{X,p)F{X,p). 



(27) 
(28) 



(29) 
(30) 



^ = (l,p/m) for non-relativistic particles and tt^ = for relativistic bosons. 
Within the same approximation level there are two alternative equations for F and F 



2n^ = V 



1 



MF - ReG^Fi, = - {{T, F} - A}) 

MF - ReG^Fo^t = \ ({r, f} - {T,^,, A}] 

with the "mass" function 

Pq — /2m — ReS'^(X, p) for non-relativistic particles. 



M(X,p) 



-w? + p^ — ReS'^(X, p) for relativistic bosons 



(31) 

(32) 

(33) 



These two equations result from sums of the corresponding Dyson's equations with their 
adjoint ones. Eqs. (^TJ) and ( P^j) can be called the mass-shell equations, since in the 
quasiparticle limit they provide the on- mass-shell condition M = 0. Appropriate combi- 
nations of the two sets (p7D-(p8|) and (pl])-(p2D provide us with retarded Green's function 
equations, which are simultaneously solved O, E3] by 



G 



R 



M(X,p) + ir(X,p)/2 



A 



ReG^ 



M2 + rV4' 

M 



(34) 



With the solution (0) for G^ equations (^7|) and ( ^8]) become identical to each other 
as well as Eqs. (p^) and (| 



I). However, Eqs. ( 27 )-( 28) still are not identical to Eqs. (31)- 
(^), while they were identical before the gradient expansion. Indeed, one can show 
that Eqs. 



( p8| ) differ from Eqs. (pll)-([3^) in second order gradient terms. This 
is acceptable within the gradient approximation, however, the still remaining difference 
in the second-order terms is inconvenient from the practical point of view. Following 
Botermans and Malfliet |8[, we express 1;^ = Tf + 0{dx) and T,^ = T{1 T f) + 0{dx) 
from the l.h.s. of mass-shell Eqs. ( pi] ) and (|52|), substitute them into the Poisson bracketed 
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terms of Eqs. (p7[ ) and ([28|), and neglect all the resulting second-order gradient terms. The 
so obtained quantum four-phase-space kinetic equations for F = f A and F = (1 =1= f)A 
then read 

P(/A)-{r/,ReG^} = C, (35) 
V ((1 T f)A) - {r(l T /), ReG^} = tC. (36) 

These quantum four-phase-space kinetic equations, which are identical to each other in 
view of the retarded relation (|33), are at the same time completely identical to the corre- 



spondingly substituted mass-shell Eqs. (^) and (|3^ . 

The validity of the gradient approximation [^] relies on the overall smallness of the 
collision term C = {gain — loss} rather than on the smallness of the damping width F. 
Indeed, while fluctuations and correlations are governed by time scales given by F, the 
Kadanoff-Baym equations describe the behavior of the ensemble mean of the occupation 
in phase-space F{X,p). It implies that F{X,p) varies on space-time scales determined by 
C. In cases where F is not small enough by itself, the system has to be sufficiently close to 
equilibrium in order to provide a valid gradient approximation through the smallness of the 
collision term C. Both the Kadanoff-Baym (KB), eq. (P?]), and the Botermans-Malfiiet 
(BM) choice ( p5D are, of course, equivalent within the validity range of the first-order 
gradient approximation. Frequently, however, such equations are used beyond the limits 
of their validity as ad-hoc equations, and then the different versions may lead to different 
results. So far we have no physical condition to prefer one of the choices. The procedure, 
where in all Poisson brackets the Fin and Fout terms have consistently been replaced by 
/F and (1 =F /)F, respectively, is therefore optional. However, in doing so we gained some 
advantages. Beside the fact that quantum four-phase-space kinetic equation ( |35[ ) and 
the mass-shell equation are then exactly equivalent to each other, this set of equations 
has a particular features with respect to the definition of a non-equilibrium entropy fiow 
in connection with the formulation of an exact H-theorem in certain cases. If we omit 
these substitutions, both these features would become approximate with deviations at the 
second-order gradient level. A numerical scheme of the BM choice in application to heavy 
ion collisions is constructed in refs. |^, |50[| . 

The equations so far presented, mostly with the Kadanoff-Baym choice (0), were the 
starting point for many derivations of extended Boltzmann and generalized kinetic equa- 
tions, ever since these equations have been formulated in 1962. Most of those derivations 
use the equal-time reduction by integrating the four-phase-space equations over energy 
Po, thus reducing the description to three-phase-space information, cf. refs. [pll, E2|, Rl, 



54, 55, 56, 57, pq, p9[ and refs. therein. This can only consistently be done in the limit 



of small width F employing some kind of quasi-particle ansatz for the spectral function 
A{X,p). Particular attention has been payed to the treatment of the time-derivative parts 
in the Poisson brackets, which in the four-phase-space formulation still appear time-local, 
i.e. Markovian, while they lead to retardation effects in the equal-time reduction. Gen- 
eralized quasiparticle ansatze were proposed, which essentially improve the quality and 
consistency of the approximation, providing those extra terms to the naive Boltzmann 
equation (some times called additional collision term) which are responsible for the cor- 
rect second-order virial corrections and the appropriate conservation of total energy, cf. 
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56[] and refs. therein. However, all these derivations imply some information loss 
about the differential mass spectrum due to the inherent reduction to a 3-momentum rep- 
resentation of the distribution functions by some specific ansatz. With the aim to treat 
cases as those displayed in Figs. 2 and 3, where the differential mass spectrum can be 
observed by di-lepton spectra, within a self-consistent non-equilibrium approach, one has 
to treat the differential mass information dynamically, i.e. by means of Eq. (|3^ avoiding 
any kind of quasi-particle reductions and work with the full quantum four phase-space 
kinetic Eq. (pS]). In the following we discuss the properties of this set of quantum kinetic 
equations. 



5 <I>-Derivable Approximations 



The preceding considerations have shown that one needs a transport scheme adapted to 
broad resonances. Besides the conservation laws it should comply with requirements of 
unitarity and detailed balance. A practical suggestion has been given in ref. |^ in terms 
of cross-sections. However, this picture is tied to the concept of asymptotic states and 
therefore not well suited for the general case, in particular, if more than one channel 
feeds into a broad resonance. Therefore, we suggest to revive the so-called ^-derivable 
scheme, originally proposed by Baym [|l^ on the basis of the generating functional, or 



partition sum, given by Luttinger and Ward |T9[, and later reformulated in terms of 
path-integrals 



30 



With the aim to come to a self-consistent and conserving treatment on the two-point 



function level, we generalized the ^-functional method [|l^, |T8l to the real-time contour 
(C) in ref. It was based on a decomposition of the generating functional F with 

bilocal sources into a two-particle reducible part and an auxiliary functional $ which 
compiles all two-particle-irreducible (2PI) vacuum diagrams 



iF{G',0,A} = iF°{G'°}+ / dxC\ct>,d^ct>} 



+ < 



-iE 



C 



(37) 



2PI 



±ln(l-0G'°0S) ±QGqE 



i<l>{G,0,A} 



Here CP{(t)) is the free classical Lagrangian of the classical field 0, G° and G denote the 
free and full contour Green's functions, while S is the full contour self-energy of the 
particles. Contrary to the perturbation theory, here the auxiliary functional $ is given 
by all two particle irreducible closed diagrams in terms of full propagators G, full time 
dependent classical fields (f) and bare vertices. Upper signs in Eq. (|37D relate to fermion 
quantities, whereas lower signs, to boson ones, while and nx count the number of 
self-energy insertions in the ring diagrams and the number of vertices in the diagrams of 
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$, respectively, A is the scaling factor in each vertex. The stationarity conditions 



6r{G,(f),X}/6G = 0, 6r{G,(f),X} /6(j) = (38) 

provide the set of coupled equations of motion for the classical fields (j) and Green's 
functions G (Dyson Eq.) 



0(a;) = 0'^(x) - dyG'^ix, y)J{y), (39) 
G(x, y) = G\x, y)+ f dzdz'G%x, z)J:{z, z')G{z', y), (40) 



where the superscript ^ marks the free Green's functions and classical fields. The func- 
tional ^{G, (j)} acts as the generating functional for the self-energy S and source currents 
J{x) via the functional variations 

The advantage of this formulation is that $ can be truncated at any level, thus defining 
approximation schemes with built in internal consistency with respect to conservation laws 
and thermodynamic consistency. For details we refer to |18, 19] and our previous paper 



||2 1|| . Note that $ itself is constructed in terms of "full" Green's functions, where "full" 
now takes the sense of solving self-consistently the Dyson's equation with the driving term 
derived from this approximate $ through relation (0). It means that even restricting 
ourselves to a single diagram in $, in fact, we deal with a whole sub-series of diagrams 
in terms of free propagators, and "full" takes the sense of the sum of this whole sub- 
series. Thus restricting the infinite set of diagrams for $ to either only a few of them or 
some sub-series of them defines a <l>-derivable approximation. Such approximations have 
the following distinct properties: (a) they are conserving, if $ preserves the invariances 
and symmetries of the Lagrangian for the full theory; (b) lead to a consistent dynamics, 
and (c) are thermodynamically consistent. These properties originally shown within the 
imaginary time formalism with a time-dependent external perturbation [T^, |19|] also to 
hold in the genuine non-equilibrium case formulated in the real-time field theory pT |. 



Transport equation (pSf ) weighted either with the charge e or with 4- momentum p'^, in- 
tegrated over momentum and summed over internal degrees of freedom like spin (Tr) gives 
rise to the charge or energy-momentum conservation laws, respectively, with the Noether 
4-current and Noether energy-momentum tensor defined by the following expressions 



nX) = eTr J ^^27T^F{X,p), (42) 



(2vr) 

e^"(X) = Try A^27r>''F(X,p) + g"" (^'°*(X) - ^p°*(X)) . (43) 



Here 



5X{x) 



(44) 

A=l 
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is interaction energy density, which in terms of $ is given by a functional variation with 
respect to a space-time dependent couphng strength of interaction part of the Lagrangian 
density C'''^ ^ \{x)C''^\ cf. ref. ^ 



The potential energy density £^°^ takes the form 



ReS^F + ReG^-F 
A 



(45) 



Whereas the first term in squared brackets complies with quasiparticle expectations, 
namely mean potential times density, the second term displays the role of fluctuations in 
the potential energy density. 

The conservation laws only hold, if all the self-energies are ^-derivable. In ref. p^ , 
it was shown that this implies the following consistency relations (a) for the conserved 
current 

"{ReS^,F} ^1^"^^ ^ 



iTr 



i2nY 



ReG 



A 



C 



0, 



(46) 



and (b) for the energy-momentum tensor 



Qu i^^pot _ ^int^ j JjReE^, f} - jReG^, + C 



(47) 



They are obtained after first-order gradient expansion of the corresponding exact relations. 
The contributions from the Markovian collision term C drop out in both cases, cf. Eq. 
(0) below. The first term in each of the two relations refers to the change from the free 
velocity v to the group velocity Vg in the medium. It can therefore be associated with a 
corresponding drag-flow contribution of the surrounding matter to the current or energy- 
momentum flow. The second (fluctuation) term compensates the former contribution 
and can therefore be associated with a back-flow contribution, which restores the Noether 
expressions (|4^ ) and (^3]) to be indeed the conserved quantities. In this compensation 
we see the essential role of fluctuations in the quantum kinetic description. Dropping or 
approximating this term would spoil the conservation laws. Indeed, both expressions (|42|) 
and (1^) comply with the quantum kinetic equation (|35|) , being approximate (up to the 
first-order gradient terms) integrals of it. 

Expressions (^2|) and (^31) for 4-current end energy-momentum tensor, respectively, 
as well as self-consistency relations (iB) and (^) still need a renormalization. They are 
written explicitly for the case of non-relativistic particles which number is conserved. 
This follows from the conventional way of non-relativistic renormalization for such parti- 
cles based on normal ordering. When the number of particles is not conserved (e.g., for 
phonons) or a system of relativistic particles is considered, one should replace F{X,p) 
^ (^F{X,p) =p F{X,p)^ in all above formulas in order to take proper account of zero point 
vibrations (e.g., of phonons) or of the vacuum polarization in the relativistic case. These 
symmetrized equations, derived from special (=f) combinations of the transport equa- 
tions ( 



) and (|36D , are generally ultra-violet divergent, and hence, have to be properly 
renormalized at the vacuum level. 
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6 Collision Term 



To further discuss the transport treatment we need an exphcit form of the colhsion term 
(PPP, which is provided from the $ functional in the — h matrix notation via the variation 
rules (0) as 



C{X,p) 



51$ 



6F{X,p) 



F{X,p) 



5FiX,p) 



F{X,p). 



(4J 



Here we assumed $ be transformed into the Wigner representation and all ^iG~~^ and 
iG^ to be replaced by the Wigner-densities F and F. Thus, the structure of the collision 
term can be inferred from the structure of the diagrams contributing to the functional 
To this end, in close analogy to the consideration of ref. p3[, we discuss various decom- 



positions of the ^-functional, from which the in- and out-rates are derived. For the sake 
of physical transparency, we confine our treatment to the local case, where in the Wigner 
representation all the Green's functions are taken at the same space-time coordinate X 
and all non-localities, i.e. derivative corrections, are disregarded. Derivative corrections 
give rise to memory effects in the collision term, which will be analyzed separately for the 
specific case of the triangle diagram. 

Consider a given closed diagram of at this level specified by a certain number nx 
of vertices and a certain contraction pattern. This fixes the topology of such a contour 
diagram. It leads to 2"^ different diagrams in the — h notation from the summation over 
all — h signs attached to each vertex. Any — h notation diagram of $, which contains 
vertices of either sign, can be decomposed into two pieces in such a way that each of the 
two sub-pieces contains vertices of only one type of sign^ 




F,...Fi. 



(49) 



with Fi ■ ■ ■ FmF[ ■ ■ ■ F^ linking the two amplitudes. The V*{X; pi, ...p[, ...) and Vg(X; pi, ...p[, 
amplitudes represent multi-point vertex functions of only one sign for the vertices, i.e. 
they are either entirely time ordered (— vertices) or entirely anti-time ordered (+ ver- 
tices). Here we used the fact that adjoint expressions are complex conjugate to each other. 
Each such vertex function is determined by normal Feynman diagram rules. Applying 
the matrix variation rules (^) , we find that the considered $ diagram gives the following 
contribution to the local part of the collision term (^9]) 



d'^Pi d^p'i 



(27r)4 {27Ty 



■■■R 



E 



"^To construct the decomposition, just deform a given mixed-vertex diagram of $ in such a way that 
all -I- and — vertices are placed left and respectively right from a vertical division line and then cut along 
this line. 
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X 



{F,...Fi... - F,...F[...} (27r)V (y^P^ - 



(50) 



with the partial process rates 



R{X;p,,...p[,...)= Re{V:{X;p,,...p[,...)Vp{X;pr,...p[,...)}. (51) 

(a/3)e4> 

The restriction to the real part arises, since with {a\f3) also the adjoint {(3\a) diagram 
contributes to this collision term. However these rates are not necessarily positive. In 
this point, the generalized scheme differs from the conventional Boltzmann kinetics. 
An important example of approximate $ which we extensively use below is 



1$ 





(52) 



where logarithmic factors due to the special features of the ^-diagrammatic technique 
are written out explicitly, cf. ref. ^2\. In this example we assume a system of fermions 
interacting via a two-body potential V = Vo6{x — y), and, for the sake of simplicity, 
disregard its spin structure. The $ functional of Eq. (|52D results in the following local 
collision term 



C 



loc 



d'^Pi d^p2 d^p3 



(27r)4 (27r)4 (27r)^ 
{P + Pi-P2- Ps) {F2F3FF1 - F^F^FFi 



(53) 



where d is the spin (and maybe isospin) degeneracy factor. From this example one can 
see that the positive definiteness of transition rate is not evident. 

The first-order gradient corrections to the local collision term ( pUD are called memory 
corrections. Only diagrams of third and higher order in the number of vertices give rise to 
memory effects. In particular, only the last diagram of Eq. ( ^2t) gives rise to the memory 
correction, which is calculated in ref. [p2|. 



7 Entropy 

Compared to exact descriptions, which are time reversible, reduced description schemes 
in terms of relevant degrees of freedom have access only to some limited information and 
thus normally lead to irreversibility. In the Green's function formalism presented here the 
information loss arises from the truncation of the exact Martin-Schwinger hierarchy, where 
the exact one-particle Green's function couples to the two-particle Green's functions, cf. 



refs. [15, 48 1, which in turn are coupled to the three-particle level, etc. This truncation is 



achieved by the standard Wick decomposition, where all observables are expressed through 
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one-particle propagators and tlierefore higher-order correlations are dropped. This step 
provides the Dyson's equation and the corresponding loss of information is expected to 
lead to a growth of entropy with time. 

We start with general manipulations which lead us to definition of the kinetic entropy 
flow [gg. We multiply Eq. by - ln(F/A), Eq. by (t) ln(F/A), take their sum, 



integrate it over d^p/(27r)^, and finally sum the result over internal degrees of freedom 
like spin (Tr). Then we arrive at the following relation 

d.sUX) = Tr/ 7^1n^C(X,p), (54) 



<oc = Tr/ ^^A^{X,p)a{X,p), (55) 



where 



a{X,p) = T[lTf] ln[l T /] - / In/, (56) 
AT 

A^{X,p) = —B^ (57) 



= A 



2 



27r^ - ME"^ 



dpa J dp^ 



(5^ 



cf. the corresponding drift term (proportional to 9^/ in Eq. (|35D). The zero-components 
of these functions, and B^, have a meaning of the entropy and flow spectral functions, 
respectively, and satisfy the same sum rule as A. If the considered particle is a resonance, 
like the A or p-meson resonances in hadron physics, the B^ function relates to the energy 
variations of scattering phase shift of the scattering channel coupling to the resonance in 
the discussed above virial limit. The value Sy^^ is interpreted as the local (Markovian) 
part of the entropy flow. Indeed, the s^^^ has proper thermodynamic and quasiparticle 
limits E^. However, to be sure that this is indeed the entropy flow we must prove the 



H-theorem for this quantity. 

First, let us consider the case, when memory corrections to the collision term are 
negligible. Then we can make use of the form ( ^0]) of the local collision term. Thus, we 
arrive at the relation 



Tr / —^ln-Cioc{X,p) 
J 27r 4 F ^ ^' 




(2vr) 

X -Fi...F{...}ln^:::^(27r)V {^P.-Y.P\) ■ (59) 

In case all rates R are non-negative, i.e. -R > 0, this expression is non-negative, since 
(x — y) ln(x/?/) > for any positive x and y. In particular, i? > takes place for all 
$-functionals up to two vertices. Then the divergence of sj^^. is non-negative 

dMX) > 0, (60) 
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which proves the if-theorem in this case with (^) as the non-equihbrium entropy flow. 
However, as has been mentioned above, we are unable to show that R always takes non- 
negative values for all $-functionals. 

If memory corrections are essential, the situation is even more involved. Let us consider 
this situation again at the example of the $ approximation given by Eq. (|52D. We assume 
that the fermion-fermion potential interaction is such that the corresponding transition 
rate of the corresponding local collision term (^) is always non-negative, so that the H- 
theorem takes place in the local approximation, i.e. when we keep only C^°'^. Here we will 



schematically describe calculations of ref. p2| which, to our opinion, illustrate a general 



strategy for the derivation of memory correction to the entropy, provided the if-theorem 
holds for the local part. 

Now Eq. (|5^ ) takes the form 

8 Am = Ti' / 1^ in f + Ti' / 1^ in f C"- (61) 

where s'^^^ is still the Markovian entropy flow defined by Eq. (^). Our aim here is to 
present the last term on the r.h.s. of Eq. (^) in the form of full x-derivative 



/ f = -^.<cm(^) + 5c^er.{X) (62) 



of some function s'^^^{X), which we then interpret as a non-Markovian correction to 
the entropy flow of Eq. ([55|) plus a correction (^Cmem)- For the memory induced by the 



triangle diagram of Eq.(|52D detailed calculations of ref. |]22| show that smallness of the 
(^Cmem, Originating from small space-time gradients and small deviation from equilibrium, 
allows us to neglect this term as compared with the first term in r.h.s. of Eq. (|5^. Thus, 
we obtain 

«oc + <em) ^ Tr I ^ In ^C'^ > 0, (63) 

which is the if-theorem for the non-Markovian kinetic equation under consideration with 
■^loc + ■^mem ^^s the propcr entropy flow. The r.h.s. of Eq. ( |63|) is non-negative provided 
the corresponding transition rate in the local collision term of Eq. (^) is non-negative. 

The explicit form of s'^^^ is very complicated, see ref. In equilibrium at low 

temperatures we get s^^^^ ~ In T that gives the leading correction to the standard 
Fermi-liquid entropy. This is the famous correction |^T], ^ to the specific heat of liquid 
^He. Since this correction is quite comparable (numerically) to the leading term in the 
specific heat (~ T), one may claim that liquid ^He is a liquid with quite strong memory 
effects from the point of view of kinetics. 

8 Pion-Condensate Phase Transition 

As a further example for the role of finite width effects we consider the phase transition 
dynamics into a pion-condensate. The possible formation of such a pion condensate in 
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dense nuclear matter was initially suggested by A.B. Migdal in his pioneering work |p. 
In realistic treatments of this problem applied to equilibrated isospin-symmetric nuclear 
matter at low temperatures T <C m.,^ the pion self-energy is determined by nucleon- 
nucleon-hole, A-nucleon-hole contributions corrected by nucleon-nucleon correlations, 
TTTT fluctuations and a residual interaction [^]. A recent numerical analysis within a 
variational method with realistic two- and three-nucleon interactions gave pc — 2po for 
the critical density of vr^, tt^, 7r° condensation in symmetric nuclear matter and pc — l-3po 
for 7r° condensation in neutron matter, with po being nuclear saturation density. 

In symmetric nuclear matter the pion condensate frequency vanishes while the magni- 
tude of condensate momentum p^, is approximately given by the nucleon Fermi momen- 
tum I |~ Pp. The critical behavior of the system is determined by the effective pion 
gap 

(64) 



mm 
P 



where the momentum p = Pc corresponds to the minimum of the gap at zero mean 



CM 



CM n 9 - 




Figure 5: Effective pion gap (|6^ ) versus nuclear density, see ref. 0. 

field = H, 1^. Fig. 5 illustrates the behavior of the effective pion gap oj'^ip^ as 
a function of the baryon density p. At low densities, ReS^ is small and one obviously 
has u;^ > 0. The dashed line in Fig. 5 describes the case where the vrvr fluctuations are 
artificially switched off and the phase transition turns out to be of second order. At the 
critical point of the pion condensation (p = pc) this value of cD^ with switched-off tttt 
fluctuations changes its sign. In reality the tttt fluctuations are significant in the vicinity 
of the critical point |^4|, |6^ . The corresponding contribution to the pion self-energy 
behaves as T/a;((y9^, p^) at T >| u;^((y9^,p^)/m^ |, and a'^(Pc) does not cross the zero hue at 
allQ. Rather there are two branches (solid curves in Fig. 5) with positive and respectively 

^Here we have used the uP quantity that aheady takes account of the pion mean field as explained 
below, cf. Eq. (p9|) and the definition of after it. 
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negative value for uj'^{p^) and the transition becomes of the first order. Calculations of 



',2/ 



refs I^J, |65|, ^ demonstrate that at p > pc the free energy of the state with to (Pq) > 0, 
where the pion mean field is zero, becomes larger than that of the corresponding state 
with a;^(Pc) < ^ finite mean field. Thus, aX p = pc the first-order phase transition 
to the inhomogeneous pion-condensate state occurs. At p > pc the state with cD^(Pc) > 
is meta-stable and the state with u^^Pc) < and 7^ becomes the ground state. 

Before we discuss a self-consistent scheme for a quantitative treatment of this problem 
we like to qualitatively explain how the instability towards pion condensation develops 
dynamically. To simplify the treatment we assume that the pion density is low (p,r ^ p) 
and further use the fact that the pion is much lighter than the nucleon {mT^/mN ~ 1/7). 
This allows us to consider the pion sub-system as a light admixture in a heavy baryon 
environment, neglecting the feedback of the pions onto the baryons. It provides nucleon 
Green's functions unaffected by the pion distribution. This very approximation was used 
in the first works exploring the possibility of the pion condensation in dense nuclear matter 

H, 0. We will use it for the pion retarded self-energy, thus neglecting the contribution 
from pion fluctuations (see dashed curve in Fig. 5). Within the above approximations 



the quantum kinetic equation (|35D for the pion distribution in homogeneous and equi- 
librated baryon environment becomes 

1. 
2' 



-TB^d'^U = Ti, - TU. (65) 



Here 5^ is defined in Eq. (E3) and all subscripts tt are omitted, except for the pion 
distribution function 

We like now to illustrate that the second branch in Fig. 5 with negative a;^ constructed 
under the assumption of vanishing mean field is indeed unstable and becomes stabilized 
by a finite mean field. The instability of the system can be discussed considering a 
weak perturbation 5/^ of the pion distribution /^°) = [exp{pQ/T) — which we assume 
equilibrated in the rest frame of the system. Linearizing Eq. ( |65D we find 



'-B,d^^6U + 5U = 0, (66) 



1 

2' 

with the solution 

8Ut,p) = 5 Up) exp [-2t/B,{p)] , (67) 

where for simplicity the initial fluctuation Sfo{p) of the pion distribution is assumed to 
be space-independent. Let us consider the case, where po —>■ and | p | — Pf- This four- 
momentum region, being far from the pion mass shell, is right the region, where the pion 
instability is expected in symmetric nuclear matter. Here the real part of the pion self- 
energy ReS^ is an even function of the pion energy po while the width is an odd function 
and proportional to po for pq — >• 0. Using the results of refs. p|, Q 2pQ — dReT,^/dpo — >■ 0, 
r = P{p)pq and j3{p) ~ nin for po 0, we get Bq = P{p)/u^{p) from Eq. (^) and 
therefore 

5Mt,po = 0,p) = 5fo{po = 0,p) exp [-2u\p)t//3{p)] . (68) 

The above solution shows that for uj"^ > initial fluctuations are damped, whereas they 
grow in the opposite case. Thus, the change of sign of uj'^{p^) leads to an instability of the 
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virtual pion distribution at low energies and momenta ^ p^. The solution (^) illustrates 
the important role of the width in the quantum kinetic description. If the width would 
be neglected in the quantum kinetic equation, one would fail to find the above instability. 

The growth of the pion distribution ^/^r is accompanied by a growth of the condensate 
field (Ptt. Due to the latter the increase of the virtual pion distribution slows down and 
finally stops when the mean field reaches its stationary value. Therefore, a consistent 
treatment of the problem requires the solution of the coupled system of the quantum 
kinetic equation (|35D and the mean field equation (|39|). In order to find the behavior of 
the virtual pion distribution one also has to include the mean field contribution to the 
pion self-energy. Considering only small mean fields we keep terms of the lowest-order in 
iPt^. Then T,^ acquires an additional contribution = = 0) + \cs \ 

where Agg denotes the total in-medium pion-pion interaction. Within the same order the 
mean-field equation becomes 



^At) = 0. 



(69) 



Here we have assumed the simplest structure for the condensate field y?^ = (pT,{t)exp{ip^r), 
where (pn{t) is a space-homogeneous real function which varies slowly in time. Also one 
should do the replacement cD^(Pc) ~^ ^"^{firyPc) = ^^(Pc) + -^efr I V^tt P in the above 
equations (|65D - (|^) for the pion distribution. 

The time dependence of can qualitatively be understood inspecting the two limits 
of small and large times. At short times the mean field is still small and one can neglect 
the Aeff<^^(t) term in Eq. Then the mean field 



<^.(t) =^.(0)exp -2u^{pJt/P{p 



(70) 



grows exponentially with time, just like the distribution function (68). Here (^7r(0) is an 
initial small fluctuation of the field. At later times the solution of Eq. (^) approaches 
the stationary limit (p-^ <^J^* with 



stat \ 2 



-^\Pc)/^cs{Pc)- 



(71) 



Since simultaneously uj'^{ipT,,pJ = oj'^{Pc) + Ks \ ^ 0, the change in the pion 

distribution 5/,^ will saturate. This stationary solution (ff^^ is stable as can be seen from 
linearizing Eq. (^9]) around iff'^^ 



^f*-5^oexp Au'{p,)t/Pip:, 



(72) 



since the exponent is negative. Here 6(po denotes an arbitrary initial space-homogeneous 
fluctuation. 

The physics can again be cast into a ^-derivable form, where the $ functional should 
include at least the following diagrams 
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+ X + >JlC +H'H > + ^ + X (73) 



Here bold and bold-wavy lines relate to baryon and pion Green's functions, respectively, 
while the wavy line terminated by cross denotes the pion condensate. Since in the broken 
phase the mean pion field mixes nucleon with A configurations we adopt the SU{4) for- 
mulation of the model, introduced in ref. There one deals with a unified description 



of baryons (A^ and A), based on 20 x 20 matrix Hamiltonian in the basis of 20 A-nucleon 
spin-isospin states. Thus, the solid lines symbolize a unified propagator matrix for A- 
resonance and nucleon. The mixing is provided by condensate-baryon coupling (diagram 

(g) ). Numerical symmetry factors are omitted in Eq. ([73|). 

Functional variation of $ with respect to propagators provides the corresponding self- 
energies. Diagrammatically this variation corresponds to cutting and opening the respec- 
tive propagator lines of the diagrams of $ in Eq. (|73p. Thus diagrams (a) to (d), (f) and 

(h) contribute to the pion self-energy. Diagram (a) accounts for the baryon particle-hole 
contributions to the pion self-energy. It includes NN~'^, AN~^, NA~^ and AA~^ terms. 
The subsequent series of diagrams (b) to (d) renormalizes baryon-pion vertex including 
baryon-baryon correlations in terms of the Landau-Migdal parameter g'. Diagram (f ) ac- 
counts for the pion fiuctuations. It is proportional to ~ T/u;((/?^, p^) and thus causes the 
transition to be of first order. This becomes especially important for the case of heated 
and even non-equilibrium dense matter, where the effective pion gap drops. One 
should notice that pion fluctuation contributions are also present in the particle-hole dia- 
gram (a) when opened perturbatively. Diagram (h) corresponds to pion interactions with 
the condensate which is responsible for the stabilization of the condensate solution ([7T|). 

Likewise cutting and opening the solid lines in $ determines the baryon self-energy 
which describes the feedback of the pions onto the baryonic subsystem. This feedback is 
required for the conserving and thermodynamically consistent treatment of the problem. 
Diagrams of the flrst line correspond to the modiflcation of the baryon motion by the 
multiple interaction with the pions corrected by correlations. Diagram (e) generates 
a purely local interaction contribution whereas diagram (g) with the coupling of the 
condensate to baryons leads to the mixing of N and A. 

Variation of $ with respect to the condensates (wavy line with cross) determines the 
source term J in the equation for the mean fleld (pOl) . The value Aefr entering Eq. (^) is 
generated by the the last two diagrams (h) and (i) of Eq. ([73|). 

The kinetic description ([35| ) for the particle distribution together with the equation 
of motion for the mean fleld (|39D is still insufficient for the numerical simulations of the 
dynamics of the phase transition. The reason is that the creation of seeds of the new 
phase, which initiate the growth of the mean fleld and the particle distribution, is due to 
fluctuations, cf. Eqs. ( |68D and ([70|). However, the scheme of Eqs. (^) and ( |39D provides 
no sources of stochastic fluctuations. Thus, it can only simulate the dynamics of one of the 
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phases rather than the transition between them. The required stochastic sources may be 
introduced into the transport theory in the spirit of the Boltzmann-Langevin approach 
developed in refs |68, |6^, 70 1 and the stochastic interpretation of the Kadanoff-Baym 



equations [^]. The stochastic transport approach offers an appropriate framework for 
the description of the unstable dynamics by means of a stochastic force in the mean field 
equation and a stochastic collision term in the transport equation, which both act as a 
source for a continuous branching of the dynamical trajectories. 

The above example shows that we really need the off-mass- shell kinetics to describe 
the dynamics of the pion-condensate phase transition, since the corresponding instability 
of the pion distribution function occurs far from the pion mass shell, cf . Eq. (|68|) . Besides 
the conserving property and thermodynamic consistency of the ^-derivable approximation 
it also leads us to the proper order of the phase transition. 



9 Summary and Prospects 

A number of problems arising in different dynamical systems, e.g. in heavy-ion collisions, 
require an explicit treatment of dynamical evolution of particles with finite mass-width. 
This was demonstrated for the example of bremsstrahlung from a nuclear source, where 
the soft part of the spectrum can be reproduced only provided the mass-widths of nucleons 
in the source are taken explicitly into account. In this case the mass- width arises due to 
collisional broadening of nucleons. Another examples considered concern propagation 
of broad resonances (like p-meson and A) in the medium. Decays of p-mesons are an 
important source of di-leptons radiated by excited nuclear matter. As shown, a consistent 
description of the invariant-mass spectrum of radiated di-leptons can be only achieved if 
one accounts for the in-medium modification of the p-meson width (more precisely, its 
spectral function). The principle of actio = re-actio was demonstrated on a pedagogical 
example when there is only vriVA coupling and in the limit of a dilute nuclear matter. We 
also expect a consistent description of chiral a-, vr- condensates together with fluctuations, 
as an immediate application of our results to multi-component systems. 

We have argued that the Kadanoff-Baym equation within the first-order gradient 
approximation, slightly modified to make the set of Dyson's equations exactly consis- 
tent (rather than up to the second-order gradient terms), together with algebraic equa- 
tion for the spectral function provide a proper frame for a quantum four-phase-space 
kinetic description that applies also to systems of unstable particles. The quantum four- 
momentum-space kinetic equation proves to be charge and energy-momentum conserving 
and thermodynamically consistent, provided it is based on a <l>-derivable approximation. 
The $ functional also gives rise to a very natural representation of the collision term. 
Various self-consistent approximations are known since long time which do not explic- 
itly use the ^-derivable concept like self-consistent Born and T-matrix approximations. 
The advantage the $ functional method consists in offering a regular way of constructing 
various self-consistent approximations. 

We have also addressed the question whether a closed non-equilibrium system ap- 
proaches the thermodynamic equilibrium during its evolution. We obtained a definite 
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expression for a local (Markovian) entropy flow and were able to explicitly demonstrate 
the if-theorem for some of the common choices of $ approximations. This expression 
holds beyond the quasiparticle picture and thus generalizes the well-known Boltzmann 
kinetic entropy. Memory effects in the quantum four-momentum-space kinetics were dis- 
cussed and a general strategy to deduce memory corrections to the entropy was outlined. 

At the example of pion-condensate phase transition in symmetric nuclear matter we 
demonstrated important role of the width effects in the dynamics and we formulated a self- 
consistent $ derivable scheme for the transport treatment of this problem. An interesting 
application of such self-consistent transport description is possible to dynamics of the 
phase transition of a neutron star to the pion or kaon condensate state accompanied by 
the corresponding neutrino burst. In view of the letter, another application concerns 
description of the neutrino transport in supernovas and hot neutron stars during first 
few minutes of their evolution. At an initial stage, neutrinos typically of thermal energy, 
produced outside (in the mantel) and inside the neutron-star core, are trapped within 
the regions of production. However, coherent effects in neutrino production and their 
rescattering on nucleons |^ reduce the opacity of the nuclear-medium and may allow for 
soft neutrinos to escape the core and contribute to the heating off the mantle. The extra 
energy transport may be sufficient to blow off the supernova's mantle in the framework 



of the shock-reheating mechanism [^. The description of the neutrino transport in the 
semi-transparent region should therefore be treated with the due account of mass-widths 
effects. 

Further applications concern relativistic plasmas, such as QCD and QED plasmas. The 
plasma of deconfined quarks and gluons was present in the early Universe, it may exist in 
cores of massive neutron stars, and may also be produced in laboratory in ultra-relativistic 
nucleus-nucleus collisions. All these systems need a proper treatment of particle trans- 
port. Perturbative description of soft-quanta propagation suffers of infrared divergences 
and one needs a systematic study of the particle mass-width effects in order to treat 
them, cf. ref. [^. A thermodynamic ^-derivable approximation for hot relativistic QED 
plasmas — a gas of electrons and positrons in a thermal bath of photons — was recently 
considered in ref. ||7^. Their treatment may be also applied to the high-temperature 
super-conductors and the fractional quantum Hall effect |7^, Formulated above ap- 
proach allows for a natural generalization of such a ^-derivable schemes to the dynamical 
case. 
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